Method and Apparatus for 
Predicting the Failure of a Component 



The patent claims priority pursuant to 35 U.S.C. § 1 19(e)l to provisional application 
5 60/243,880 filed 10/26/2000. 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 
This invention relates to a method and apparatus for predicting failure of a 
10 component. More specifically it relates to a method and apparatus for predicting failure of a 
component using a microstructure-based fatigue model developed using probabilistic 
methods. 

DESCRIPTION OF THE RELATED ART 
Ever increasing economic demands require that components and structures have 
15 high reliability. It is usually the case that similar components experience large variations in 
their total life. For example, one machine element may last many years, yet another element 
produced by the same manufacturer may fail in a few months. Traditional design methods 
account for large uncertainty or scatter in life by using design safety factors. Low 
probabilities of failure can usually be ensured only by applying large safety factors during 

20 

component design. But safety factors can add unnecessary redundancy, weight, and cost to 
a component because the true safe life of the component may be much greater than the 
predicted safe life. Furthermore, to ensure safety, components are often retired well before 
the useful lifetime is exhausted to insure that no component fails during operation (i.e. 
Aircraft parts). Designing cost effective and highly reliable structures and maximizing part 
^ life therefore requires the ability to accurately assess a component's safe life. 

Assessing a component's safe life often requires determining the effects of fatigue on 
a component. Fatigue can occur in any device that has moving components. Fatigue can 
also occur in cases where the movement is imperceptible as do, for example, bridge 
elements or railroad tracks. Components also often fail in an insidious manner, giving no 
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u prior indication that damage has occurred. In the case of an aircraft engine, for example, 
such fatigue failures can be catastrophic. 

Failure analysis has revealed that actual component loadings are often well below 
the steady loads that can cause failure. What distinguishes these failures is the fact that the 
loads have been applied a large number of times. Fatigue, more specifically fatigue crack 
initiation and growth, is therefore a damage mechanism that degrades the reliability and safe 
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life of components subjected to repeated loads. Designing for fatigue is difficult because it 
typically exhibits a large variation in its effect on similar components. 

In the case of fatigue failure, scatter in life is indicated by a coefficient of variation 
(COV) which is usually determined based on a wide range of fatigue life tests on many 
5 similar specimens. Even under well controlled laboratory tests of annealed smooth 
specimens at room temperature the COV varies from less than 10% to over 500% for 
different steel alloys. Thus, the considerable scatter in the fatigue reliability of components 
in operation may be substantially attributed to considerable scatter of component material 
fatigue behavior. 

10 Life scatter is due to the fact that, in general, materials have inhomogeneous 

microstructures. To the naked eye, it may appear that a material is composed of continuous 
homogeneous material. Microscopic examination reveals, however, that metals, for 
example, are comprised of a discontinuous inhomogeneous material consisting of individual 
crystalline grains, pours, and defects. Cracks nucleate and grow on the order of grain size 

15 according to the properties of the individual grains, with growth rate as varied as grain 
properties. As these cracks grow the rate and behavior of the crack approaches the bulk or 
average properties of the material. Therefore, for large cracks, traditional crack growth 
methods are appropriate. Traditional methods, however, fail to determine the probability of 
crack initiation or to describe crack growth of near-grain-sized cracks. In many applications 

or* 

failure can occur before the fatigue damage reaches the long crack stage because the energy 
associated with the damage is very high although the damage is very small. 

Current fatigue life prediction methods in metallic components consider three 
stages: crack initiation, long crack propagation, and final fracture. Long crack propagation 
and final fracture are stages of damage accumulation that are well characterized using 
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computational models such as linear elastic or elastic-plastic fracture mechanics. 

Crack initiation is the early stage of damage accumulation characterized by small 
cracks; cracks with depths less than several grain diameters. These have been observed to 
deviate significantly from predicted long crack fracture mechanics. The deviation is 
attributed to the heterogenous material in which small cracks evolve. The crack initiation 
JKJ phase accounts for the majority of scatter in fatigue life for many alloys. The crack 
initiation stage contains two phases: crack nucleation and small crack growth. 

Crack nucleation is a locally complex process of crack formation on the 
microstructural scale. One example of a crack nucleation mechanism is the smooth fracture 
at angles inclined to the loading direction that is exhibited by materials having a propensity 
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for planar slip. Crack growth is the similarly complex process that occurs once a crack has 
been nucleated. 

Current crack initation models are based on empirical testing. This causes crack 
initiation models to be simple parametric functions of applied stress variables. As such, 
^ these macrostructural models assume the material to be homogeneous and continuous. 
Statistical concepts have been used to develop empirical fatigue life models where the 
independent variable, usually applied global stress or strain, is considered deterministic and 
the dependent variable, usually life, is considered random. But these models do not account 
for the mechanisms of the microstructural parameters that regulate fatigue damage. Since 
10 this major source of scatter ( i.e. microstructural crack initiation, which includes both crack 
nucleation and small crack growth) is not included in these models, they are necessarily 
unsatisfactory because they cannot represent the heterogeneous material in which the 
damage processes occur. 

Because traditional crack initiation models are empirical, they cannot represent 
15 conditions not included in the established database test program. Sequential variation is one 
such condition. Sequential variation is due to the component being used for many different 
f J scenarios. For example, an automobile is driven differently during each trip and an 

unlimited number of sequential variations would have to be considered. While current long 
crack growth propagation models are able to account for sequential variation in the 

20 

component usage, it is not practical from a time and cost standpoint to include sequential 
variations in the applied stress under most test programs. Therefore crack initiation tests are 
generally conducted at a maximum i.e., "worst case", stress to ensure safety. Thus 
components are generally systematically over-designed assuming unlikely worst case 
material properties to compensate for the lack of true understanding of the material fatigue 
behavior. Additionally, these traditional models are unsatisfactory for predicting individual 
component failure because the vast majority of components by definition do not possess 
these unlikely worst case material properties. Predictions made using such models are 
based upon the worst case material properties and are thus inaccurate. 

Figure 1 depicts three levels of fatigue damage that may occur in a typical high 
strength component. First, a crack nucleates 200 on a small scale on the order of the grain 
size. Then the crack grows as a microscopically small crack 202 in which the crack lies in 
relatively few grains. The material properties, averaged along the front of the crack, 
approach bulk or average material properties as the crack grows and the number of grains 
intercepted by the crack front increase. When the material properties intercept enough 
grains traditional crack growth techniques such as linear elastic fracture mechanics 204 may 
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be applied. But the majority of crack life is spent in the nucleation and small crack growth 
regime for high strength components. Thus, understanding the early crack behavior is most 
important. As a result, there exists a need for a method and apparatus for accurately 
predicting component failure that account for the microstructural properties of materials and 
sequential variation in the loading, and relate them to fatigue scatter. 



BRIEF SUMMARY OF THE INVENTION 
This invention provides a method for predicting the failure of a component. A 
preferred embodiment of the method comprises obtaining a Finite Element Model (FEM) of 
) the component; analyzing the FEM to obtain stresses at its nodes; identifying a subset of the 
nodes as significant nodes based on the stresses; determining a Representative Volume 
Element (RVE) for the significant nodes; developing an RVE microstructure-based failure 
model for the RVEs; simulating a component life using RVE microstructure-based failure 
models to produce a result related to the component life; performing the simulating a 
* plurality of times to produce results related to the component life; preparing statistics using 
the results; and comparing the statistics to a probability of failure (POF) criteria to 
determine whether the performing predicted failure for the component. 

The invention also provides an apparatus for predicting the failure of a component. 
A preferred embodiment of the method comprises a central processing unit (CPU); an 
1 output device for displaying simulated fatigue results; an input device for receiving input; 
and a memory comprising: instructions for receiving input comprising a component's 
material characteristics; instructions for using RVE microstructure-based failure models and 
the input and predicting failure of the component, the predicting comprising: simulating a 
component using at least one RVE microstructure-based failure model, the simulating 
producing a result related to component life; performing the simulating a plurality of times 
to produce results related to component life; preparing statistics using the results; and 
comparing the statistics to a probability of failure (POF) criteria to determine whether the 
performing predicted failure for the component; and instructions for displaying a result from 
the predicting. 

In addition the invention provides a method for determining the orientation factor for 
a grain slip system of a material comprising: obtaining equations that relate a stress 
direction to a material's potential slip systems; simulating a grain orientation of the material, 
the simulating comprising: using probabilistic methods to generate a slip plane normal 
angle for each of the potential slip systems; inputting the normal angle into the equations to 
obtain a potential orientation factor for each of the potential slip systems; and selecting the 
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least of the potential orientation factors as a grain orientation factor for the simulated grain 
orientation; repeating the simulating for a defined number of grains and obtaining a plurality 
of grain orientation factors; and creating a statistical distribution of the plurality of grain 
orientation factors to determine an orientation factor for the grain slip system. 

Finally the invention provides an apparatus for determining the orientation factor for 
a grain slip system of a material comprising: a central processing unit (CPU); an output 
device for displaying simulated fatigue results; an input device for receiving input; and a 
memory comprising: instructions for receiving input; instructions for simulating a grain 
orientation of the material, the simulating comprising: relating a stress direction to a 

10 material's potential slip systems with equations; using probabilistic methods to generate a 
slip plane normal angle for each of the potential slip systems; inputting the normal angle 
into the equations to obtain a potential orientation factor for each of the potential slip 
systems; and selecting the least of the potential orientation factors as a grain orientation 
factor for the simulated grain orientation; instructions for repeating the simulating for a 

15 defined number of grains and obtaining a plurality of grain orientation factors; and 

instructions for creating a statistical distribution of the plurality of grain orientation factors 
to determine an orientation factor for the grain slip system. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

20 

The foregoing and other aspects and advantages of the present invention will be 
better understood from the following detailed description of preferred embodiments of the 
invention with reference to the drawings, in which: 

Figure 1 illustrates three levels of fatigue damage in a high strength component; 

Figure 2 is a microscopic view of a metallic structure showing grains of various 
2 ^ shapes and sizes; 

Figures 3(a) - 3(e) depict a flowchart of a preferred embodiment of a method of the 
invention; 

Figure 4 is a diagram of an apparatus for predicting the failure of a component that 
incorporates a preferred embodiment of the present invention; 

30 

Figure 5 depicts a flowchart of preferred embodiment of a method of the invention 
for determining a grain orientation factor; 

Figure 6 is a diagram of an apparatus for determining the orientation factor for a 
grain slip system that incorporates a preferred embodiment of the present invention; 

Figure 7 depicts a crack tip slip band in multiple grains in Example 1; 
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Figure 8 depicts a random microstructure generated by a Monte Carlo simulation in 
Example 1; 

Figure 9 is a drawing of a wheel with five holes for Examples 2 and 5; 
Figure 10 is a close-up of a Finite Element Model showing a region of high stress 
5 and the location of the FEM nodes in that high stress region for Examples 2, 3, and 5; 

Figure 1 1 is a graph of the S-N curve for smooth round bar specimens of the wheel 
material for Example 2; 

Figure 12 is a close-up of a Finite Element Model near a high stress region for 
Example 2; 

10 Figure 13 is a graph showing the statistical distribution of fatigue predicted by the 

probabilistic microstructure-based failure model for the hole of Example 2; 

Figure 14 is a histogram of fatigue lives for 1000 small size bars tested at 90 ksi for 
Example 3; 

Figure 15 is a graph showing the distribution of medium size bars tested at 90 ksi for 
15 Example 3; 

Figure 16 is a graph showing the relationship between the nodal von Mises stress 
and the normalized hole diameter for Example 4; 

Figure 17 is a graph showing the probability density function of the fatigue life of 
the failed holes for the deterministic and random holes of Example 4; 

20 

Figure 18 is a graph showing the probability of failure (POF) cumulative density 
function of the fatigue life for wheels without correlation in Example 5; 

Figure 19 is a graph showing the probability of failure (POF) cumulative distribution 
function of the fatigue life for wheels with various amount of correlation in Example 5; 
Figure 20 illustrates the orientation of a slip system in Example 6; 

25 

Figure 21 depicts the stereographic triangle of a uniform statistical distribution of the 
Schmid factors for randomly oriented grains generated by a Monte Carlo simulation in 
Example 6; 

Figure 22 depicts the stereographic triangle of the statistical distribution of Schmid 
factors for texture oriented grains in Example 6; 

30 

Figure 23 is a top view of the slip plane of Figure 19 in Example 6; 
Figure 24 is a side view of the slip plane of Figure 19 in Example 6; and 
Figure 25 is a graph illustrating the predicted probability density function of the 
reciprocal Schmid factor for slip in fully lamellar titanium aluminide (TiAl). 

Like reference numerals refer to corresponding elements throughout the several 
35 drawings. 
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DETAILED DESCRIPTION OF THE INVENTION 
The present invention uses a virtual prototyping technique that relies on computer 
simulation of real material behavior to predict when a component will fail due to fatigue. 
The fatigue life of a component is recognized to be a characteristic property of the material 

5 of which the component is composed. A component is made up of an ensemble of discrete 
microscopic structural elements of a material such as grains, colonies, and nodules. Many 
different materials are composed of such discrete microscopic structural elements, such as 
metals, intermetallics, ceramics, and some plastics. The method of the invention is aimed at 
predicting the variation in fatigue life based on the statistical variation of the microscopic 

3 structural of the material. Material parameters at metallic grain level are used along with 
fundamental physics-based models to predict the damage as it accumulates from the 
nucleation of cracks, through small crack growth and long crack growth, to final failure. 

The computer simulates many "identical" components but uses a different sample of 
material microstructure for each simulation. The microscopic structure of each simulated 

) material model or "realization" for each component is properly sampled from the known or 
specified range of material microstructures. Each of the elements is then virtually tested 
using computer simulation to simulate real-world usage conditions. The virtual testing 
allows data to be produced quickly on thousands or even millions of components. This 
virtual testing addresses variation in the microscopic substructure, illustrated in Figure 2, by 

' modeling the grain size, grain orientation, micro-applied stress and micro-yield strength as 
random. These parameters are then used in modeling crack nucleation and small crack 
growth. All of the variation in the long crack growth is simulated by the variability in crack 
growth rate coefficient. 

Preferred embodiments of the present invention use probabilistic analysis methods 
to address the effects of the random nature associated with material microstructure, loading, 
and manufacturing differences, assuming that the underlying physics of the component 
behavior are deterministic and that the random nature of the component response is 
attributed to the variability in the input to the component and the parameters defining the 
failure physics. 

Now referring to Figures 3(a) - 3(e), which contain a flowchart of a preferred 
embodiment of a method of the invention, the method generally includes identifying the 
material microstructure 1, identifying the damage interaction with the material 
microstructure 2, and developing a microstructure-based failure model 3 (Figure 3b). 

Referring to Figure 3(a), identifying the material microstructure 1 includes: 
conducting literature search 4; making direct observations about a materials microstructural 
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arrangement based on specimen samples 5; and evaluating microstructure geometry 
statistics 6 based on the literature search and direct measurements of specimens. These 
statistics will vary depending on the material and microstructural arrangement but this will 
usually include grain size, orientation, and volume fraction estimates. Identifying the 

5 material microstructure 1 culminates with characterizing the physical grain microstructure 7. 
Identification of the damage interaction with the material microstructure 2 includes 
processes necessary to characterize active damage mechanisms 15 that become the basis of 
microstructure-based failure model development 3 (Figure 3(b)). First a literature review 8 
is undertaken to determine if information exists about either the bulk elastic material 

3 characteristics 14 or the mechanical (i.e. loaded) microstructural characteristics 9. Defect 
properties assessment 10 then defines grain slip planes, pores, or inclusions which are likely 
to cause local plastic deformation. There are a wide variety of other potential defects within 
any given material. Crack nucleation properties 1 1 defines the mechanisms that cause local 
plastic deformation to nucleate cracks. Short crack growth properties 12 then defines the 

> active mechanisms at the short crack tip that govern the erratic behavior of short cracks, 
such as grain boundary blockage, grain orientation, and the local frictional strength. Long 
crack growth properties 13 defines crack growth rate parameters as well as threshold 
characteristics. Determining bulk elastic material characteristics 14 encompasses a number 
of properties including: shear modulus, Poisson's ratio, and specific fracture energy; 

) although the appropriate properties vary greatly depending upon the material and 
microstructural arrangement. Then, after determining bulk elastic material characteristics 
14 and mechanical microstructure characteristics 9 it is possible to characterize the damage 
mechanisms 15. 

Now referring to Figure 3(b), development of the microstructure-based failure model 
1 3 is generally accomplished by determining the stages of damage accumulation 16, 
developing crack nucleation models 22, developing short crack growth models 30 (Figure 
3(c)), developing long crack growth models 38 (Figure 3(c)), and linking (i.e. sequencing 
and/or nesting) the crack nucleation, short crack growth, and long crack growth models to 
produce an overall failure model 46 (Figure 3(c)). 

Determining the stages for damage accumulation 16 begins by determining the 
number of simultaneous damage mechanisms 17 that exist (i.e. ductile nucleation in one 
phase occurring simultaneously with brittle nucleation in a second phase). Then a 
mechanism sequencing or nesting strategy 18 is developed that, in general, links these 
models. Linking is needed to predict the fatigue response at the component level because 
fatigue damage can start at a very small scale and grow to final fracture. The lower level 
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model uses the appropriate parameters to determine the initial state of the next level. The 
next level uses the results from the previous level along with the appropriate parameters 
specific to its level to determine the initial state of the next level and so on. Through the use 
of these "nested" modeling techniques, component reliability can be linked to the physical 
5 mechanisms at each level. Additionally, by modeling each level of the fatigue process 
individually, and rigorously linking the levels, various size effects are included. Next make 
an evaluation as to the number of nucleation models 19, short crack growth models 20, and 
long crack growth models 21 required. Each of the mechanisms characterized 15 (Figure 
3(b)) will have a corresponding model. Multiple models may be needed for each level due 
10 to multiple mechanisms for each damage level, and different mechanisms may be active 
simultaneously 17. 

Next develop the crack nucleation models 22 determined appropriate for the 
mechanisms of damage accumulation 19. First, identify the important random variables 23 
based on the required model inputs 24, or outputs 25 from other models (should the 
15 nucleation model be nested). Then relate the important variables to crack nucleation 26 
through literature review 27 and development through theoretical considerations 28. Crack 
life and crack cycles are similar concepts when discussed regarding fatigue failures. Finally, 
define model output 29 (parameters that could also serve as inputs to other models at the 
next level). 

20 

Now referring to Figure 3(c), develop the short crack growth models 30 determined 
appropriate 20. The method is analogous to developing the crack nucleation models 22. 
Identify important random variables 31 based on required model inputs 32 or outputs from 
other models 33 (in the case where the short crack growth model is nested). Then relate the 
important variables to short crack growth 34 through literature review 35 and development 
25 through theoretical considerations 36. Finally, define model output 37 (parameters that 
could also serve as inputs to other models at the next level). 

Similarly, develop the long crack growth models 38 determined appropriate 21. 

Identify important random variables 39 based on required model inputs 40 or outputs from 

other models 41 (in the case where the long crack growth model is nested). Then relate the 
30 • 

important variables to long crack growth 42 through literature review 43 and development 
through theoretical considerations 44. Finally, define model output 45. 

After developing the appropriate models, these nucleation, short crack, and long 
crack growth models are sequenced or nested 46 following the strategy developed 18 
(Figure 3(b)) to produce the overall microstructure-based failure model. 

35 
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Microstructure-based failure model 46 is limited to a single microstructure and 
single loading conditions. Most real-world components will have many material 
microstructural arrangements and, more likely than not, experience multiple loading 
conditions. Now referring to Figure 3(d), to apply the method to a real-world component, 
5 develop or obtain a conventional Finite Element Model (FEM) 48. Then analyze the FEM 
and obtain the stress at each node 50. With these stresses, identify the significant nodes 52. 
Significant stresses could be those above material fatigue strength or to a safety factor 
applied to fatigue strength as determined by one of ordinary skill in the art. Around these 
significant nodes develop a Representative Volume Element (RVE) 54. An RVE is a finite 
0 region of the FEM that has a consistent stress. The information for each RVE should 
include the stress on the RVE, the dimensions, and the properties and microstructure of the 
material within the RVE. Determine the microstructures within each RVE and apply the 
process of Figures 3(a) - 3(c) to each RVE 56 to obtain an overall microstructure-based 
^ failure model 46 (Figure 3(c)) for each RVE. At some point prior to beginning the first 
3 simulation of the component, identify a probability of failure (POF) convergence criteria 58 
for the eventual statistics that describe the simulated component fatigue failures 80 (Figure 
3(e)). The necessary statistical significance will be different for different situations. 

The component is simulated by using the overall microstructure-based failure model 
developed for each RVE 56. First establish the density of potential nucleation sites within 
) each RVE using the probabilistic-based Monte Carlo (MC) method 62. Then establish the 
number of potential nucleation sites within each RVE using MC 64. The use of MC 
methods for establishing the number and density of potential nucleation sites is documented 
within the open literature. 

Next determine the cycles to failure for each potential nucleation site 66. First use 
' MC methods to establish values for each of the random variables 68 in the overall failure 
model 46 (Figure 3(c)) for the potential nucleation site in question. The appropriate set of 
those values is input into the microstructure-based failure model for the potential nucleation 
site in question 70 to determine the cycles to failure for that site 72. The process is repeated 
for each potential nucleation site within the RVE 72 and the "life" of the RVE is established 
74. The life of the RVE is the smallest number of cycles to failure for any of the included 
potential nucleation sites. 

Now referring to Figure 3(e), repeat simulating the component 60 (Figure 3(d)) for 
each RVE 76 in this simulated component. The predicted life of this component 78 is the 
minimum number of cycles to failure for any RVE. Include this result with any previous 
predictions and determine the statistics on these component fatigue failures 80. For 
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example the statistics could include a POF, a cumulative distribution function (CDF), or a 
probability density function (PDF) of failure characteristics. One of ordinary skill in the art 
will recognize that the desired component statistics will vary considerably based on the 
problem at hand. Next, check the statistics against the POF convergence criteria 82 and, if 
5 not met, return to simulate another component failure 60 and repeat the process. 

Should the RVE be so large that computer analysis of every potential nucleation site 
would be cumbersome, a modification can be made to the process. Step 66 would be 
modified to evaluate only a statistically significant number of potential nucleation sites and 
probabilistic methods (also called system reliability methods) would be used to estimate 
10 RVE life 74, rather than directly computing the shortest life of each potential nucleation 
site. 

Another aspect of the invention is that it can be modified to account for an FEM that 
accounts for variation in stress and/or dimension. Here, step 50 is modified so the FEM 
analysis results in a statistical distribution of stresses rather than a single value of stress. 
15 This statistical distribution of stresses may be found experimentally, or developed using any 
number of probabilistic methods. Then in step 68, stress will also be one of the random 
variables whose value is established using MC analysis, also a probabilistic method. 

Still another modification of the invention adds a spatial correlation to the FEM. 
This correlation is beneficial when a component has multiple locations of similar geometric 
detail. An optional step finds the fatigue life for regions of the component with similar 
geometric detail 75, repeats the process for each RVE 76 and then determines the 
component fatigue life based in part on a spatial correlation from the information gained 
from step 75. The spatial correlation comes from using a common established life for the 
RVE 74 whenever encountering an RVE that is one of the group having similar geometric 
25 detail. This may lead to the use of the various probabilistic methods to calculate component 
life (based on the particular circumstances). An RVE may, in fact, only be two dimensional, 
but it is nevertheless referred to as a volume element. 

The probabilistic method used in determining cycle to failure for each nucleation site 
66, providing values for random variables 62 (Figure 3(d)), 64, 68, and estimating RVE 
u fatigue life 74 include Fast Probability Methods (FPM) and Simulation Methods (ST). FPM 
techniques include response surface FPM and direct FPM. Direct FPM methods will always 
provide a solution, but when a response surface may be used its use can increase the 
efficiency of the prediction calculations. A response surface, however, cannot be formed 
when considering variables that vary with time and, thus, present discontinuities. Direct 
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FPM are then necessary, although such variables may possibly be handled using multiple 
nested response surface equations, a single response surface equation will not suffice. 

Other potential FPM approaches include First Order Reliability Methods (FORM), 
Second Order Reliability Methods (SORM), Advanced Mean Value (AMV) methods, and 
5 Mean Value (MV) methods. Potential ST approaches include Monte Carlo (MC) and 
importance sampling. MC methods are used in this preferred embodiment of the invention 
for simulating components using the microstructure-based failure models for the RVEs 60. 

Now referring to Figure 4, which is a diagram of an apparatus for predicting the 
failure of a component that incorporates a preferred embodiment of the present invention, 

10 apparatus 91 generally includes an input device 90, a central processing unit (CPU) 92, an 
output device 94, a bus 98, and a memory 96. Memory 96 contains instructions for 
receiving input 100, simulating a component life 102, preparing statistics 112, and 
displaying the prediction 116. Input 100 contains the component's material characteristics 
and values for other variables necessary to predict failure, such as the number of 

15 components to simulate, etc. Simulating a component life 102 entails using a micro- 
structure-based failure model to simulate the life of RVEs by establishing the density and 
number of nucleation sites within an RVE 104, determining the cycles to failure for each 
nucleation sites within the RVE 106, establishing an RVE life 108 for each RVE based on 
the cycles to failure for each nucleation site, and predicting the life of the simulated 

20 

component 1 10 based on the established RVE lives. Once that prediction is made for that 
single simulated component, that prediction is added to the group of any previous 
component life predictions and statistics prepared 1 12 to describe the group of predictions. 
These statistics are compared 114 to a previously input Probability of Failure (POF) criteria, 
and if not met another component life is simulated 102. If the POF criteria are met then the 

25 prediction is displayed 1 16. 

Another cause of the large scatter in the fatigue response is the limited ductility 
caused by the ordered microstructure that constrains plastic slip deformation to a limited 
number of crystallographic directions. Describing the slip behavior is essential to 
developing computational material models that predict mechanical properties such as 

30 fracture, fatigue, and creep behavior. The present invention provides a method for 
statistically describing the grain orientation factor for slip. Deformation is strongly 
orientation dependent for many high strength alloys. With the method of the present 
invention, the orientation dependence of the deformation can be described statistically and 
probabilistic microstructure-based models can be used to predict the scatter in the 
macrostructural material response such as fatigue. 
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Now referring to Figure 5, which contains a flowchart of preferred embodiment of a 
method of the invention for determining a grain orientation factor, the first step is to identify 
the active slip systems 140. Then the geometry between relative slip systems is determined 
142 to obtain equations that relate direction and potential of the slip systems 144. Next the 
5 material texture is determined 146 and the number of grains to be simulated defined 148. If 
a material is untextured it implies that it exhibits a uniformly random grain orientation. 
Simulating grain orientation 150 includes randomly generating an angle of slip plane normal 
using MC 152, inputting the generated angle of slip plane normal into equations (from step 
144) 154 to generate an orientation factor ("M") for each potential slip system 156. The 
^ minimum orientation factor generated is selected as the grain orientation factor for this 
simulated grain 158. Simulating grain orientation 150 is repeated for the user-defined 
number of grains 160 to obtain a grain orientation factor for each simulated grain. These 
grain orientation factors are statistically analyzed to predict the grain orientation factor for 
the material 162. 

15 Now referring to Figure 6, which is a diagram of an apparatus for determining the 

orientation factor for a grain slip system that incorporates a preferred embodiment of the 
present invention, apparatus 93 generally includes an input device 90, a central processing 
unit (CPU) 92, an output device 94, a bus 98, and a memory 96. Memory 96 generally 
contains instructions for receiving input 118, predicting a grain orientation factor for the slip 

20 system 119, and displaying the prediction 130 of the distribution of the grain orientation 
factors for the slip system. Predicting a grain orientation factor for the slip system 119 in 
turn includes simulating a grain orientation 120 of a material in question by using 
probabilistic methods to generate a slip plane normal angle for each potential slip system, 
inputting the normal angle into equations 122 that relate the directions of the potential slip 

^ systems to obtain a potential orientation factor for each potential slip system, selecting the 
minimum potential orientation factors 124 as the grain orientation factors for that simulated 
grain orientation, repeating the simulation 126 of the grain orientation for a defined number 
of grains and obtaining a grain orientation factor from each simulation, and creating a 
statistical distribution 128 of the grain orientation factors to determine the orientation factor 

30 for the grain slip system. 

The principals of the present invention are further illustrated by the following 
examples. These examples describe possible preferred embodiments for illustrative 
purposes only, they do not limit the scope of the invention as set forth in the appended 
claims. 

35 
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EXAMPLE 1 

The following example applies the method of the present invention to a y prime- 
strengthened nickel superalloy to derive a microstructure-based fatigue failure model and is 
illustrated by Figures 7 and 8. It also reflects the numbering and methods of Figures 3(a) - 
3(e) throughout. The example is applicable to the general class of nickel superalloys. The 
discrete microscopic structural elements of a material were identified 1 through literature 
search 4 and direct microscopic observations 5. Through literature search 4 it was 
determined that the alloy was face center cubic single phase with polycrystalline grains. For 
this example, no direct observations of specimen samples 5 were required. Information in 
the literature was used to estimate the statistics defining the geometry of the microstructure 
6 as characterized 7 by the values in Table 1 for the grain diameter, slip orientation of 
surface grains, and slip orientation of interior grains. 



Table 1: Characterization of the physical microstructure. 



Variable 


\ Description 


Distribution 
Type 


Distribution 
Parameters 


Average 


cov 


d 


Grain diameter 


Lognormal 


Ap=-0.076 £=0.39 


0.005m 


0.40 


M s 


Schmid 
Orientation Factor 


Curve Fit 


2.21 


0.08 


M T 


Taylor Orientation 
Factor 


Curve Fit 


3.07 


0.13 



Identifying the interaction of the fatigue damage with the material microstructure 2 
started with literature review 8. Dislocation moving on slip planes and pinning at the grain 

25 boundaries served as the initial defects 10 to accommodate plastic deformation. Crack 
nucleation 11 was determined to occur when the accumulated dislocations exceeded the 
materials specific fracture energy. The crack may not grow, i.e. the crack growth may 
arrest, depending of the size of the crack and the applied load. 

The short cracks were determined to grow 12 by emitting dislocations from the crack 

30 tips along the slip planes of the grains ahead of the short crack. The dislocation movement 
caused a zone of plastic deformation at the crack tip. The zones propagated freely if the 
crack tip was far from a grain boundary. As the crack tip approached a grain boundary the 
zone pinned at the grain boundary. The crack may not grow into the next grains depending 
on the size of the plastic zone. The size of the plastic zone was dependent on the crack size, 

35 the applied load, and the size, orientation and strength of the grains surrounding the plastic 
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zone. The crack was considered a short crack until the crack is of sufficient size such that 
the crack tip intercepted sufficient grains to assume bulk properties. Long crack growth 13 
was described using conventional methods based on linear elastic fracture mechanics as is 
typical in the open literature. Review of the literature 8 identified suitable bulk elastic 
^ material properties 14 from empirical microstructural studies reporting bulk properties. 
The microstructural fatigue failure model 3 was determined based on the 
observations of the damage interaction with the microstructure 2. Because the example 
material was single phase and the only defects were dislocations accumulating on naturally 
occurring slip planes, no mechanisms acted simultaneously 17. The linking strategy 18 was 
10 that the crack nucleated within a grain, then the crack grew as a short crack until it reached 
sufficient size, then the crack grew as a long crack. A single model was used for crack 
nucleation 19. Two models were used for short crack growth 20. One short crack model 
was used when the slip band tip was free to propagate and the other when the slip band tip 
m was blocked by a grain boundary. A single model was used for long crack growth 21 . 
C3 i:> The crack nucleation model 22 was used to determine the number of cycles needed 

Sy to crack the grain. The crack nucleation size was set to be equal to the grain size. The 
W important variables 23 governing the crack nucleation were: W { the plastic energy cause by 
|p the i th fatigue cycle; G, the shear modulus; t, the local applied stress; k, the factional stress 
! which must be overcome to move dislocations; v, Poisson's ratio; d, the grain diameter; a, 

|2 the local applied normal stress; W s , the specific fracture energy per unit area; and M, the 
1 1 orientation factor. For this example, all the important variables were required inputs 24 
Q there were no outputs 25 from other models. 

Li. 

r " An example of a crack nucleation model that related the important variable to the 

crack nucleation cycles and size 26 for a planar slip driven metallic grain was developed 
25 from theoretical considerations 28 as well as literature review 27 and was described 
mathematically in Equation 1: 



Eq. (1) 

(Act 



30 

Wi = 



-2k 



M _ 
AG 



n{\ - v)d 



35 
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The plastic energy was found for each cycle and total plastic energy for each grain 
was the summation of the plastic energy from each cycle and described by Equation 2: 



Eq. (2) 

5 Q 

Wt = ; 
i=i 

where W T is the accumulated fracture energy for Q cycles. The number of cycles needed to 
nucleate a crack the size of the grain was determined to be when W T exceeded W s . The 
10 variables, which may have had significant local variation, were: W s , a, M, k, and d. The 
output 29 from the crack nucleation model was N n , the number of cycles needed to grow a 
crack to the size of the grain. 

The short crack growth model 30 was used to determine the number of cycles 
needed to grow the crack from the initiation size until traditional long crack growth 
15 mechanics 38 could be assumed. The important variables 31 governing the crack nucleation 
were: cp„the crack tip opening displacement (CTOD); and C',an experimentally determined 
coefficient. The CTOD was a function of the applied load and the microstructural 
»j parameters surrounding the crack tip such as: G, the shear modulus; k 9 the factional stress 
t which must be overcome to move dislocations; n, Poisson's ratio; d, the grain diameter; a, 

20 the local applied normal stress; W s9 the specific fracture energy per unit area; and M, the 
f orientation factor. These important variables were required inputs 32. One other important 
i variable, an output from other models 33, was a n , the crack size at the end of the crack 
nucleation period. 

An example of a short crack growth model that related the important variable to the 
25 short crack growth cycles and size 34 for a planar slip driven metallic grain was developed 
from theoretical considerations 36 as well as literature review 35 and was described 
mathematically as follows in Equation 3, assuming a simple linear relationship between the 
short crack growth rate and the CTOD: 



30 Eq. (3) 



da „ /a 
— =C A<p, 

dN 



where: a was the crack size; N was the number of cycles; C" was an experimentally 
determined coefficient; and q>, was the CTOD. The CTOD was predicted based on the 

35 
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interaction of the crack tip plastic zone with the microstructure. The crack size at any 
number of cycles Q was defined 37 by Equation 4: 



Eq. (4) 



Q 



a Q = a n + ^da 



where a n was the crack size at the end of the crack nucleation period. 

An example of a CTOD model for a planar slip driven metallic grain was described 
10 mathematically as follows. Considering a crack with the crack tip in the f grain and the 
slip band tip in the n th grain as shown in Figure 7. If the slip band was propagating (not 
blocked by the grain boundary), the size of the slip band zone could be found from Equation 
5: 



where: x, was the applied resolved shear stress in the i' h grain; k, was the frictional stress of 
the i th grain; a was the crack length; c was the crack plus slip band length; and L, was the 
distance from the free surface to the grain boundary of the i th grain preceding the slip band 
tip as shown in Figure 7. 
25 The CTOD was defined by Equation 6: 



15 Eq- (5) 



0 = 




arcsin — 
c 




Eq. (6) 



2k a 




30 



35 




- 17- 



CAl -288004.1 



A = G/2k (I - v) for edge dislocations 
A-G/2x for screw dislocations 



5 For the blocked slip band, the size of the slip band zone was defined by Equation 7, the 
CTOD defined by Equation 8, and the microscopic stress intensity factor at the slip band tip 
by Equation 9. 



Eq.(7) 

10 

Eq. (8) 



(0= L-a 



nA n l A a 



15 



^cL^) 



2k ] a 
p = 1 arccos — 



7ZT } C 



20 A fe-Uj^ 

" " arccos 



J 



4-1 



V c 



Eq. (9) 

25 



To account for the crystallographic orientation of the individual grains, the applied 
resolved shear stress x in Eq. 5 through 9 was replaced with: 

(7 

30 T = — — for surface gains 
M s 

a 

1 - "77— for interior gains 

7* 
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where a was the local normal applied stress, M s was the reciprocal Schmid factor and M T 
was the Taylor factor. 

Predicting the CTOD using Eq. 6 or 8 required determining the microstructural 
properties along the crack path. Considering a random array of grains as shown in Figure 8, 
5 these properties were determined to be as follows. A crack nucleated in the surface grain X 0 
and then grew as a semi-circle through zones in which the effective material properties were 
uniform. The boundaries of the zones were represented by the concentric half circles. The 
zones were composed of grains represented by the semi-circular segments. The arc length 
of the semi-circular segments was a random variable equal to the grain diameter. The 
10 surface grains were represented by the intersection of the zones and the surface. 

After successful crack nucleation, the crack grew from grain X 0 into zone 1. In the 
example shown in Figure 8, zone 1 contains three grains. The surface length 11 of zone 1 is 
the simple arithmetic average of the grain diameters shown by Equation 10. 

15 Eq. (10) 

d u +d l2 +d l3 
3 

The effective material property Pleff of zone 1 is the average of the properties of 
20 the individual grains P u weighted with the area of the grain. P leff represents the local 
factional strength k or the local applied stress, a, and is shown as Equation 11. 



Eq. (11) 



25 



P = 



dn+d? 2 +d? 2 



30 



In the n th zone composed of j grains, the surface length is given by Equation 12 and 
the effective material property by Equation 13. As the crack becomes long, I approaches the 
mean grain size and P eff approaches the bulk properties. 



Eq. (12) 



35 
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Eq. (13) 

p _ P * idni 

5 

Using the concepts of effective material properties, crack growth was modeled as 
one-dimensional. Considering a cut along the x-axis (Section A-A in Figure 8), the fatigue 
damage was modeled as a one-dimensional crack growing through zones of varying size /„ 
and varying effective material properties P eff . 

1 0 After determining the material properties along the crack path the CTOD was 

determined using Eq.6 or 8. Eq. 3 was used to determine the incremental crack growth with 
each application of loading N. The output of the short crack growth model 37 was the short 
crack growth life (the number of cycles for the crack to growth through the microstructure 
of Figure 8) and the short crack size Z f . 

1 J The long crack growth model 38 was used to determine the number of cycles needed 

to grow the crack from the short crack size until crack reached the critical crack size. The 
important variables 39 governing the crack nucleation were: C, the Paris coefficient; n, the 
Paris exponent; Ac, the globally applied stress range; and p, the shape factor (1.12 4n for a 

20 half penny shaped crack). These important variables were required inputs 40. Another 

important variable, an output from other models 41, was a sc , the crack size at the end of the 

short crack growth period. 

An example of a long crack growth model for a planar slip driven metallic grains 

was described mathematically using a model from the literature 43. In this case no 
25 additional theoretical considerations 44 were required. The incremental crack growth rate in 

the long crack growth stage was calculated using Paris law and is given by Equation 14, 

with the crack size at any number of cycles given by Equation 15. 



Eq. (14) 



30 



Eq. (15) 



35 
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The output 45 from the long crack growth model was the Ng, the long crack cycle to 
reach the critical crack size. 

The overall failure model 46 was produced by summing number of cycles for the 
three phases namely, crack nucleation, short crack growth, and long crack growth. 

5 

EXAMPLE 2 

The following example applies a preferred embodiment of the method of the present 
invention to predicting fatigue failure of a bolt hole using probabilistic microstructure-based 
models using Finite Element Method (FEM) results. It reflects the numbering of Figures 

10 3(a) - 3(e) throughout and is further illustrated using Figures 9 - 13. This example 

illustrates the development of methods to model geometrically complex components with 
complex stress distributions and predict the failure of the component. The output from 
typical FEM analyses are used to find the elastic global stress at each node 50 of the FEM 
model. The macro-stress does not consider any microstructural interactions as opposed to 

15 the microstress that is the stress caused by the interaction of the microstructure. A 

representative volume (or surface area) element (RVE) is determined for each node 54. The 
RVE grid is different from the FEM grid. Actually, not all nodes of the FEM model need be 
considered. The RVE development is performed for each node with a stress above the 
fatigue limit 52. In most cases only a few of the nodes need to be analyzed. 

^® Each RVE is then modeled as a simple test specimen using the method of Figures 

3(a) to 3(c) to provide a probabilistic microstructure-based failure model 46 of the RVE. 
Each element can have a different microstructure, loading, failure mechanisms, surface area 
and volume. The component is modeled as an assemblage of RVEs 60. Each RVE has a 
stress and size and microstructural based failure model. 

Figure 9 shows a wheel with 5 holes. The wheel has a constant thickness of 0.5 
inches. The holes have a diameter of 0.08 inches. The wheel transmits 2700 pound of force 
to a shaft bolted to the holes. A finite element model is developed 48 of the loaded wheel 
and the von Mises stress predictions 50 are determined. The loading condition causes high 
stress around the holes. Figure 10 shows a close up of the high stress area along with the 

30 location of the nodes of the FEM model for one hole. The stress level is represented by 
different shades in Figure 10. The stress associated with each shade is shown at the right of 
the figure in units of pounds per square inch (psi). The maximum stress is 51215 psi at 
node 124. 

35 Table 2: Input to Example 2. 
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Variable 


Description 


Distribution 
Type 


Distribution 
Parameters 


Average 


COV 


C 


Paris Law 

V^UClllClCIll 


Lognormal 


0.034 


C=0.30 


4.0X10^ /Win 


0.30 


C 


CTOD Law 
Coefficient 


Deterministic 


N/A 


0.10 


N/A 


d 


Grain diameter 


Lognormal 


0.076 


C=0.39 


.005w 


0.40 


' *comp 


# of component 


Deterministic 


N/A 


100,00 


N/A 


G 


Bulk Shear 
Modulus 


Deterministic 


N/A 


nxioV 


N/A 


k 


Frictional Strength 


Weibull 


n=l.l 

2 


Q=3.7 


I0,000psi 


0.30 


r M 

K cHt 


Critical 
Microstructural 
Stress Intensity 

Factor 


Deterministic 


N/A 




N/A 


M s 


Schimd 
Orientation Factor 


Curve Fit 


2.21 


0.08 


M T 


Taylor Orientation 
Factor 


Curve Fit 


3.07 


0.13 


N 


Paris Law 
Exponent 


Deterministic 


N/A 


3 


N/A 


W s 


Specific Fracture 
Energy 


Deterministic 


N/A 


3600 lbs/in 


N/A 


V 


Poisson's Ratio 


Deterministic 


N/A 


0.3 


N/A 


a 


Applied Micro- 
stress 


Normal 




a=0.3 


15000 psi* 


0.30 


*Note: The stress level is also the applied global stress. 



30 To perform a typical traditional fatigue analysis, assume the S-N curve of Figure 1 1 

is available for the wheel material. This curve was developed from smooth round bar (SRB) 
test based on the model 46 and the material input properties of Tables 1 and 2. The SRB's 
had a surface area of 0.61 in 2 . The "min" line in Figure 1 1 refers to the -3a probability of 
failure. The fatigue limit is 40 ksi. A typical fatigue life analysis finds the allowable cycles 

35 
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from the S-N curve for the maximum predicted stress from the FEM model. Using Figure 
11, a stress of 51215 psi yields a median life of 17600 cycles and a -3a life of 2500 cycles. 

A probabilistic microstructural fatigue analysis uses much of the same information 
as the typical fatigue analysis. What is not included in typical fatigue analyses is the size 
5 effect. The size effect consideration is critical to reliability analysis. In reality, the typical 
fatigue analysis above predicts the fatigue life of a SRB stressed to 51.2 ksi. Similitude is 
assumed i.e., it is assumed that the wheel is similar to the test specimen. Similitude can 
only be achieved if the test specimen and the wheel have a similar geometry and stressed 
area. Although the size effect phenomenon is generally accepted and near component shape 
10 specimens are sometimes tested, often only standard test specimen data is available. 

Probabilistic microstructural fatigue models address the size effect issue by using the 
primitive physics that governs the underlying fatigue behavior 2. By using primitive physics 
models, probabilistic microstructural models predict the test specimen behavior. Typical 
fatigue models are curve fits to the test specimen behavior. Their bases lies in the test 
15 specimens fatigue behavior. Thus, the assumption is made that the wheel is similar to the 
test specimen. This assumption is incorrect but it is still used to approximate the fatigue 
response of the wheel. 

Although probabilistic microstructural models use specimen test data to identify 
fundamental physical phenomenon 5, their basis lies at the microstructural level. The 

20 

similitude assumption made with the probabilistic microstructural models is that the 
microstructure of the wheel is similar to the microstructure of the test specimens. 
Microstructural similitude is often easier and less costly to achieve than similitude at the 
specimen level. 

To apply the probabilistic microstructural fatigue model, the area at each stresses 
level must be determined 54. Since the nodal stresses (or strains) are most often used in 
fatigue analysis (as opposed to the average element stresses), the area near each node is 
assumed to have consistent stress equal to the nodal stress. The area calculation is 
performed for each significant node with a stress above the fatigue limit 52. The fatigue 
limit is assumed to be 40 ksi. Five nodes had a stress level above the fatigue limit. Figure 

30 

12 shows a close-up of FEM model near the high stressed region. Element numbers are at 
the center of the elements. Node numbers are at the corners of the elements. Consider node 
124 that lies in element 3, 4 and 47. The FEM analysis calculates the area of each element. 
Element 3 had 3 nodes. One third of the area of element 3 is attributed to node 124. 
Likewise for element 4. Element 47 had 4 nodes. One fourth of the area of element 47 is 
attributed to node 124. The areas attributed to node 124 at each element are then summed. 
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This is the area of a RVE 54 that is used in the probabilistic fatigue model. Table 3 shows 
the stress and area for each RVE of the five node above the fatigue limit for the hole. 

In this example assume that each RVE consist of statistically similar microstructure. 
The probabilistic microstructural fatigue model 46 and the material input properties of 
5 Tables 1 and 2 are applied to each RVE 56. It is determined that 10,000 Monte Carlo 
simulations are needed for statistical significance 58. Each grain is a potential nucleation 
site. The probabilistic microstructural fatigue model 46 includes the information on the 
grain size 62. The grain size is random. Monte Carlo simulation is used to choose the 
number of grains needed to cover the RVE area 64. One grain is chosen at random along 

10 with the properties associated with the grain 68. 

The microstructural fatigue analysis 70 is performed to determine the number of 
cycles needed for the crack to grow from this grain to failure 72. The process 68 through 72 
is repeated for each grain of the RVE to determine the cycles to failure for each nucleation 
site 66. The fatigue life of each of the RVE is determined 74 as the minimum life 

15 associated with each grain. This process 62 through 74 is repeated for each of the five RVEs 
76. 

The life of the hole 78 is equal to the life of the RVE with the minimum life. The 
location of the RVE with the minimum life is also recorded. This is the location of the 
initiation site. Ten thousand holes are simulated in this manner 82. 

20 

Table 3: RVE of the hole model. 



Node 


Stress (ksi) 


Area (sq in) 


Initiation site (%) 


117 


45.0 


1 2.53E-02 


7 


118 


41.4 


1.20E-02 


0 


120 


47.5 


1.79E-02 


13 


122 


46.1 


1.71E-02 


7 


124 


51.2 


1.37E-02 


36 



The histogram and cumulative density function 80 of the fatigue life for the hole are 
shown in Figure 13. The median fatigue life is 386,000 cycles. The minimum allowable 
fatigue life (minus three standard deviation life) is 6000 cycles. The percent of holes that 
initiated failure at each location is shown in Table 3. 

This is more than the traditional fatigue analysis life of 2500 cycles because the 
highly stressed area of material around the hole was smaller than the 0.61 in 2 that was the 
underlying assumption in the traditional analysis. The median life is much higher than the 
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traditional fatigue analysis life of 17,600 cycles because 37% of the holes did not fail at the 
test suspension (10 10 cycles). Traditional fatigue analysis cannot predict test suspensions. 
Either all of the components fail or none of the components failure in a traditional fatigue 
analysis. Traditional fatigue analysis also cannot differentiate initiation sites. 

5 

EXAMPLE 3 

The following example illustrates a preferred embodiment of the method of the 
invention that uses probabilistic microstructure-based models using FEM results to large 
RVEs. The example is further illustrated by Figures 10, 14, and 15. The probabilistic 

10 microstructure fatigue model evaluates the fatigue resistance of each nucleation site of the 
RVE. A large structure such as a bridge could easily have very large RVEs that could have 
more than 10 10 nucleation sites. If, during the fatigue simulation, we wish to consider 
thousands of structures, the number of nucleation sites to consider becomes large and the 
computational time can become significant. In this example, we will extend the concept of 

15 the extreme value statistics to fatigue of large RVEs. 

Consider a large RVE that has a fatigue response controlled by transgranular (non- 
defect) crack nucleation. The fatigue response of the structure is an extreme value problem 
because the grain that nucleates a crack and grows to failure first will control the fatigue 
life. The grain that caused failure is the extreme grain. 

20 

Extreme value statistics can be extended to surface area (or volume if subsurface 
nucleation is an active failure mode) by determining the fatigue life of a small area and then 
consider the large RVE as being compose of many small areas. This means that in Figure 
3(d) the life of the RVE 74 is not based on determining the cycles to failure for each 
nucleation site 66 but simulating a statistically significant sample of nucleation sites and 

25 

approximating the life of the RVE 74 based on the extreme value statistics of the sample. 

In the example of the bridge RVE, assume the RVE has an area of 10 in 2 . The 
statistical distribution of fatigue life for a small area of 1 in 2 is found. The life of the 10 in 2 
RVE is the minimum life of 10 of the small area samples. Using area as the basis for our 
extreme value distribution has an additional benefit. The fatigue failures may be nucleated 

30 

at defects. The defect distributions are usually defined by area or volume. In the discussion 
below, extreme value statistics is applied at the grain level and then extended to surface 
area. 

Considering that there are n grains on the surface of the test specimen. The life of 
the specimen is given by Equation (16). 

35 
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Eq.(16) 



Y is the minimum life andX 7 s are the fatigue lives of the individual grains in the specimen 
72. It is assumed that the fatigue of the individual grains, X l9 . . . ,X n are independent and 
5 identically distributed. Based on the above assumption, the probability distribution of Y can 
be derived from the initial distribution of Xj where X i? is the statistically significant sample 
of nucleation sites (grains). 

It is observed that if F, the smallest among (X 1? . . ,XJ is greater than some value, y, 
then all other values in the sample must be greater than y also. Hence, the probability that Y 
10 is greater than y can be written as Equation (17). 

Eq. (17) P(F > y)=P(X ] > y, X 2 > y,..., X n > yf 

Equation (17) can be rewritten as: 

15 

Eq.(18) UFy(y) = [1-Fx(y)] n 

where, FyQ is the cumulative distribution function (CDF) of 7 and F X Q is the CDF of X. 
Rearranging the terms in Equation (18), we get the CDF of Y as: 

20 

Eq. (19) F Y (y) = 1 - [1 - [1 - F x (y)Y 

The corresponding probability density function (PDF) of Y can then be determined from 
Equation (19) as Equation (20). 

25 

Eq. (20) f Y (y) = n- [1 -F&^ffy) 

where,£() is the probability density function (PDF) of X. Equations (19) and (20) are the 
exact forms of the extreme value distribution of the minimum life Y. 

30 

The above concept determines the probability distribution of the minimum life from 
a sample of n grains. This can simply be extended to determine the probability distribution 
of the minimum life from a sample of n surface areas. For example, X t could represent the 
minimum life of a representative surface-area-element of say 0.01 in 2 . Then for a bar with a 
surface area of 1 in 2 , n (the total number of representative surface-area-elements in the bar 
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or the sample size) would be 100. Equation (20) is used to determine the minimum fatigue 
life for the bar 78. 

As an example, consider the alloy modeled using 46 and the material input 
properties of Tables 1 and 2. Let us simulate the fatigue test of 1000 area of 1 in 2 bars of 
this alloy and a load of 90 ksi. Figure 14 shows the histogram of the simulated test results. 

We can use Eq. (19) and the simulated fatigue distribution from the small (1 in 2 ) 
area to estimate the fatigue distribution of a RVE with an area of 10 in 2 as shown in Table 4 
where the first column is the number of cycles for which the probability is to be determined, 
the second column is the percent of small areas that failed and the third column is the 
predicted percentage of failures for the RVE. 



Table 4: Prediction of probability of failure for RVE tested at 90 ksi. 



20 



! Cycles 

y 


% failed 
0.001 in 2 bar 

FJy) 


Predicted % failed | 
0.01 in 2 bar 
F y (y)=l-[l-F x (y)r 


50000 


0.2% 


1.98% 


60000 


1.0% 


9.56% 


70000 


2.3% 


20.7% 


80000 


3.6% 


30.6% 


90000 


4.8% 


38.1% 









To evaluate the accuracy of the predicted results in Table 4, a fatigue test of 10,000 
RVEs was simulated. Figure 15 compares the results of the simulated fatigue life 
distribution of RVE (designated as Actual in Figure 15) with the fatigue life distribution 
predicted using Eq. (19) and the results of the simulated fatigue life distribution of 1000 
areas of 1 in 2 (designated as Predicted). In general, the predicted results are in close 
agreement with the actual results. 

Determining the statistical distribution of the 10.0 in 2 RVE base on the 1 in 2 area 
specimen require about one tenth the computation time as direct simulation of the 10.0 in 2 
RVE. Once the statistical distribution of the RVE fatigue life has been determined as in 
Figure 15, Monte Carlo simulation is used to establish the life for an individual RVE 74. 

EXAMPLE 4 
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The following applies a preferred embodiment of the method of the present 
invention that uses probabilistic microstructure-based models using probabilistic FEM 
results. The example is further illustrated by Figures 16 and 17. The prior wheel examples, 
Examples 1, 2, and 3, applied a deterministic FEM model. But the approach is general and 

5 valid for probabilistic FEM models as well. 

It is assumed that the wheel has random geometry and boundary conditions. The 
random variables are the wheel thickness T 9 the force F, and the diameter D of the holes. 
The deterministic FEM model described previously can be considered as a realization of the 
random model with each random variable defined at their mean value as shown in Table 3. 

) The deterministic model is extended to a random model by allowing variation in each of the 
random variables with coefficients of variation show in Table 5. For this example problem, 
all variables are assumed to be normally distributed, but the algorithm developed is not 
restricted to normal distributions. 



Table 5: Random variable statistics 



20 



Random Variable 


Mean 


Coefficient of Variation 


Thickness (T) 


0.5 in. 


0.05 


Load(L) 


2700 lbs 


0.10 


Diameter(£>) 


0.8 in. 


0.04 



Variations in the random variables cause variation in the nodal stresses of the FEM 
model. This means that the stress obtained at each node 50 will not be a described by a 
single deterministic value but by a statistical distribution. 

25 To apply the probabilistic microstructural fatigue model to a probabilistic FEM 

model, a response surface equation is determined for each RVE that relates the nodal 
stresses to the random variables. A virtual design of experiment is used to determine the 
response surface relationship. It is assumed that the random variables are physically 
independent which is reasonable for this example. A series of seven finite element analysis 

30 are made. One analysis with each of the three random variables at the mean value and two 
"perturbed" analysis for each random variable. In three of the perturbed analyses, one 
random variable is increased slightly above the mean value while the other random variables 
remain at the mean value. For the other three perturbed analyses, one random variable is 
decreased slightly below the mean value while the other random variables remain at the 

35 mean value. 
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Because a linear FEM analysis was performed 50, a direct linear relationship exists 
between the nodal stress and the applied load (F). A direct linear relationship also exists 
between the nodal stress and the wheel thickness (7). Because the relationship between the 
^ nodal stress and F and T is known, the perturbed FEM analysis for F and 7" need not be 
5 performed. The only "unknown" relationship is between the nodal stress and the hole 
diameter (D). 

Figure 16 shows the results of the hole diameter perturbed FEM models for each of 
the nodes in the highly stresses area. (The hole diameter has been normalized by dividing 
by the mean value.). 

The least square method was used to fit a second order response surface equation to 
the data in Figure 16 for each node. Knowing the second ordered relationship between the 
nodal stress and the hole diameter and the linear relationship between the nodal stress and 
the other random variables, a response surface of the form of Equation (21): 

; Eq.(21) a = tf (Ad 2 + Bd +C) 

is developed where o is the nodal von Mises stress, t is the normalized thickness f 77 T ) 
/is the normalized force (F/F mean ) and d is the normalized hole diameter. A, B, and C are 
the coefficients of the second ordered relationship between 0 and d. Table 3 shows the 
stress relationship along with the nodal area for each node. 



Table 6: RVE of the hole model. 



25 



Node 


Stress (ksi) 


Area (sq in) 


117 


tf (-12160^ -45024d +102154) 


2.53E-02 


118 


tf(48672d 2 -149700d +142378) 


1.20E-02 


120 


#"(60096^ -121096d + 108502) 


1.79E-02 


122 


tf(-62240d 2 + 5060rf + 103236) 


1.71E-02 


124 


//X-367488d 2 + 630072d - 21 1369) 


1.37E-02 



The probabilistic microstructural fatigue model 46 is applied to the wheel bolt hole 
by simulating a hole consisting of RVEs for each node in Table 3. ( The material 
microstructural fatigue property of each of the RVEs is described by the model 46 and the 
material input properties of Tables 1 and 2.) The parameters for each hole are determined by 
generating a random value for T, L, and D using the statistical distributions in Table 5. The 
nodal stresses are then determined for each RVE using the equations in Table 6. The fatigue 
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life of each of the RVEs is determined 74. The life of the hole is equal to the minimum life 
of each of the RVEs 78. Monte Carlo simulation is used to simulate many holes in this 
manner. 

A fatigue test of 10,000 holes with random T, L, and D was simulated and compared 
^ with simulated test of 10,000 holes with deterministic (nominal) macrostructural parameters 
(T, L, and D\ Forty two percent of the random holes did not fail at the test suspension (10 10 
cycles). Thirty seven percent of the deterministic holes did not fail at the test suspension. 
The probability density functions of the fatigue life for the failed holes are shown in Figure 
17. 

^ The median fatigue life of the failed deterministic holes is 80,000 cycles. The 

minimum allowable fatigue life (minus three standard deviation life) is 5300 cycles. The 
median fatigue life of the failed random holes is 35,000 cycles. The minimum allowable 
fatigue life is 2800 cycles. For the example problem, the variation of the FEM model 
parameters significantly lowers the minimum allowable fatigue life, 
y ^ The analysis presented in this example showed that probabilistic microstructural 

r|| models can be integrated with stochastic finite element models. In the example problem, 
W component geometry and applied loading were allowed to vary. The variation in the FEM 
p model parameters culminated in a variation of the nodal stresses. Response surface methods 
I , were used to determine the relationship between the nodal stress and the random FEM 
|g 20 moc jel variables. The variation in the nodal stress served as an input to the probabilistic 
?|| model 50. The hole was modeled in the probabilistic microstructural fatigue model 46 by 

Q assembling the appropriate RVEs. 

N ; 

EXAMPLE 5 

25 The following example illustrate a preferred embodiment of the method of the 

invention that uses probabilistic microstructure-based models using probabilistic FEM 
results with correlated random variables. It is further illustrated using Figures 9, 10, 18, 19, 
and 20. This example illustrates the treatment of the statistical correlation that may exist 
between similar geometric details within the same component. 

DKJ Most parameters in a structure are distributed in space, not concentrated at a point. 

Examples of such parameters are distributed loads, geometric properties such as thickness 
of a wheel, etc. Such quantities cannot be expressed as single random variables, but only as 
a collection of many random variables, which are statistically correlated. Statistical 
correlation between random variables is generally expressed by a correlation coefficient p. 

35 The value of this coefficient can vary between 0 and 1. A correlation coefficient value of 0 
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indicates that there is no correlation between the two variables and they are statistically 
independent of each other. Conversely, a correlation coefficient value of 1 indicates that 
there is perfect correlation between the random variables. 

Consider many wheels each with five boltholes as in Figure 9. The diameter of the 
5 boltholes for one wheel may be correlated because all of the holes on the wheel may have 
been drilled using the same operator, drill bit and drill press. Another wheel could have 
been machined using a different operator, bit and press. This would cause more similarity 
between diameters of the holes of a single wheel than between the holes of different wheels. 
Therefore, the diameters of the holes of a wheel are correlated. Similarly, because each 
> wheel is forged independently, the variation of wheel thickness may be correlated. The 
loads at the boltholes of a wheel are likely to be highly correlated because the entire wheel 
experiences the same service duty 

For the example problem (see Figure 9), wheel thickness, hole diameters and loads 
at the five holes are considered as random variables. The value of these random variables 
' are correlated at each hole of the wheel i.e., if the thickness at one hole is thin, the thickness 
at the other holes on the wheel will tend to be thin. The correlation among the wheel 
thickness at the five hole locations is represented by the correlation matrix shown in Table 
7. The matrix is symmetric. The matrix indicates that the correlation coefficient between 
the thickness values at any two locations is 0.7. This implies that if the thickness at any 
hole location is high, it is 70% likely that the thickness values at the other 4 hole locations 
will be high also and vice-versa. Table 7 shows that each random variable is perfectly 
correlated with itself. By using the full correlation matrix of Table 7, the correlation 
between any two locations can be defined (i.e., because hole 1 is close to hole 2, the 
correlation between 1 and 2 may be higher than the correlation between 1 and 3). If the 
correlation between 1 and 3 were 0.5, and the same relationship existed for the other holes, 
the correlation matrix would be as in Table 8. 



Table 7: Correlation matrix for uniform correlation. 
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t2 


t3 


t4 


t5 


tl 
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0.7 


0.7 


0.7 


0.7 


t2 
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0.7 


0.7 


0.7 


t3 






1 


0.7 


0.7 


t4 


sym 


met 


ric 
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0.7 


t5 










1 
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Table 8: Correlation matrix for nonuniform correlation. 



5 





tl 


t2 


t3 


t4 


t5 


tl 


1 


0.7 


0.5 


0.5 


0.7 


t2 




1 


0.7 


0.5 


0.5 


t3 






1 


0.7 


0.5 


t4 


sym 


met 


ric 


1 


0.7 


t5 










1 



A Monte Carlo simulation algorithm is used for generating correlated normal 
variables for input 24, 32, and 40 into the probabilistic microstructure fatigue program. For 
the correlation matrices shown above, the algorithm would be applied to generate sets of 5 
correlated variables. 

1 . Five uniform random numbers are generated between 0 and 1 and stored in the vector 
{u}. 

2. The vector uniform of random numbers is transformed to a vector of uncorrected 
standard normal variables, {S}. The transformation is performed using the Inverse CDF 
Method (Ang and Tang, 1984). Standard normal variables have a Gaussian distribution 
with a mean value of 0 and a standard deviation of 1. 

3. The vector {S} is transformed to vector of correlated standard normal variables, {5"} as 
follows: 

Eq.(22) {S>}= [L]{S] 



where [L] is the lower triangular matrix obtained from the Cholesky decomposition of the 
correlation matrix (Tables 7 or 8) of the 5 variables. 

4. The correlated standard normals, {5"}, are transformed to correlated normal thickness 
variables, ft J, as follows: 



Eq.(23) {*}={$'} 

30 

where flsj is the diagonal matrix of the standard deviation of the variables and {u} is the 
vector of the mean values of the variables. 

The above steps are used to determine the correlated values of the wheel thickness T, 
the applied force F and the diameter D to be used in the response surface Eq. (21) (of 
Example 4) as follows. 

35 
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Assume that the FEM model of the hole developed in the previous examples is a 
nominal hole in the wheel with five holes. Figure 10 shows a close up of the high stress 
area along with the location of the nodes of the FEM model. As in the previous examples, 
the hole is assumed to be an assembly of five RVEs described by the nodal stresses and 
5 nodal areas associated with nodes 117, 1 18, 120, 122, and 124. These were the only nodes 
that have stresses above a pre-defined fatigue limit 52 of 40 ksi. 

Assuming the wheel has random geometry and boundary conditions, the random 
variables are the wheel thickness T, the applied force F and the diameter D of the hole. The 
random variable statistics are shown in Table 3 but for this example, the random variables 
10 are correlated. 

To apply the probabilistic microstructural fatigue model to a stochastic FEM model, 
a response surface equation (21) was determined for each node that relates the nodal stresses 
50 to the random variables 56. Table 6 shows the stress relationship along with the nodal 
area for each node. 

1 The fatigue life of a hole is simulated by finding the fatigue life of each RVE 60 and 

setting the fatigue life of the hole equal to the minimum fatigue life of all of the RVEs. 
Likewise, the fatigue life of a wheel with multiple holes is the minimum fatigue life of each 
of five holes. 

A fatigue test was simulated for 10,000 wheels with random t, I, and d at each hole 

20 

and no correlation between the random variables. (This means that all off diagonal terms of 
the correlation matrix have a value of zero.) The probability of failure (POF) cumulative 
density function (CDF) of the fatigue life is shown in Figure 18. If there is no correlation 
between the random variables describing the holes of a wheel with five holes, the POF of 
the wheel can be determined from Equation (24) 

25 

Eq. (24) POF, = 1 - (1 - POF A ) 5 

POF w is the POF of the wheel and POF A is the POF of one hole. The development of Eq. 
(24) was presented in Equations (16) through (19), but Equation (24) is not valid if the 

30 

random variables are correlated. 

A simulation of 10,000 wheels was performed with the random variables fully 
correlated (p = 1.0) at each hole. Figure 19 shows that the POF of the wheel is less for 
correlated variables than for uncorrelated variables. Two more simulations were performed, 
one with p = 0.7 the other with p = 0.5. Figure 19 shows that the POF increases as p 
DO decreases. 
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EXAMPLE 6 

The following example illustrates a preferred embodiment of the method of the 
present invention for statistically describing the grain orientation factor for slip, depicted in 
> Figure 5, by applying that method to polycrystalline gamma titanium aluminide (TiAl). It is 
further illustrated using Figures 20 through 25. The geometric description of the active slip 
systems 140 for a single grain (colony) was used to determined a minimum Schmid factor 
for the grain. A Monte Carlo simulation was then used to generate an ensemble of grains 
with a uniform random distribution. 

' Since plastic flow occurs on crystallographic slip planes in particular directions, the 

local shear stress ,x, is a function of grain orientation and the applied stress. Consider a 
single grain with slip on a single plane in a single direction along x as shown in Figure 20. 
The resolved shear stress is given by Equations (25). 

Eqs. (25) T = ms 

m = cos C, cos X 

s is the applied stress, C, is the angle from the slip plane normal to the loading axis, and X is 
the angle from the slip direction to the loading axis, m is referred to as the Schmid factor, 
x is a maximum at ^ = X = 45 °, m = 1/2. x is zero when the tensile axis is perpendicular to 
the slip plane (X = 90°) or parallel to the slip plane (£ = 90°). The orientation factor is 
conventionally described by the reciprocal Schmid factor, M as Equation (26). 



Eq. (26) 



m cos^cosA 

In general, casting operations allow grains to solidify with uniformly random 
orientations. This orientation can be considered untextured. The orientations of the grains 
of an untextured (as-cast) polycrystalline material can be expressed as a uniform distribution 
of points within the stereographic triangle as shown in Figure 21. The angle of the slip band 
normal £ is uniformly distributed between 0 and 90 degrees. 

Deformation processes such as forging and rolling tend to alter the orientation of the 
grains. This preferred orientation is called textured microstructure. Micrographs of the 
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etched surface along with process simulation models can be used to determine grain 
orientation after deformation operations. The texture 146 can be described as a weighted 
probability function within the stereographic triangle as shown in Figure 22. The 
orientation with the most likely occurrence has the highest value. The orientation with the 
least likely occurrence has the smallest value. The statistical distribution of the slip factor 
can be determined using the described process except the uniform distribution between 0° 
and 90° of the slip plane normal for an untextured microstructure is replaced by the 
weighted probability function between 0° and 90°. This allows the determination of the 
effect of texturing 146 on the orientation factor for any polycrystal. 

The crystallographic orientation of gamma titanium aluminide microstructure consist 
of alternating layers of plate-like domains. There are domains of various phases within the 
grain. The domain that is most important for fatigue type loading can be considered similar 
to face center cubic (fee). The important slip system of this domain has the slip plane and 
the slip direction parallel to the plate. This system can be represented by the single slip 
system of Figure 20. The angle between the load direction and the slip plane normal £ can 
be any angle between 0° and 90°. But the minimum value of X is the complement of £ such 



The titanium aluminide grains are composed of colonies with many lamellae domains. 
There is a matrix-twin relationship between the lamellae. Each colony contains three 
matrix-twin pairs for a total of six easy slip directions, all 60° degrees apart in the same slip 
plane. Because slip occurs in both directions, there are 3 independent directions for slip; 0°, 
60° and 120°. Therefore, once one domain's orientation is defined, there will be other 
domains in the same colony with relative slip directions of 60° and 120° from the first slip 
direction, all in the same slip plane. This matrix-twin relationship effects the statistical 
distribution of the reciprocal Schmid factor. 

Consider Figure 23 where the ellipse represents the top view of the slip plane shown in 
Figure 20. Let us define an angle a of a slip direction in the slip plane such that when a = 0, 
X = 90° and when a = 90°, X = 90° - y. Figure 24 represents the side view of this 
configuration with the slip plane perpendicular to the plane of the paper and a = 0° into the 



that: 



Eq. (27) 



20 




35 
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paper. If a, represents the angle of the slip direction in the slip plane of an arbitrary slip 
system, then within the same colony two other slip systems exist a 2 and a 3 such that: 



Eqs. (28) 



= oil + 60° 
a, = a, + 120 c 



The angles are bounded such that: 



10 



Eqs. (29) 



0° <: a, < 90° 
60° < a 2 < 150° 
120° <: 03 < 210° 



r\ 

l y 
W 

CP 

o 



in 

£ 



A relationship between the angles exist such that: 

90° 



15 Eq. (30) 



<Xi = 



(90° -At) 



where X, as the angle between the slip direction and the loading axis corresponding to a v If 
we know X 1 and £ of one slip system we can use Equations (27) and (26) to find ^ and ^ for 
the other two slip systems available within the same colony such that: 



20 



25 



Eqs. (31) 



A 2 = 90°-^H. 
90° 



/U = 90°- 



180 o -a 2 
90° 



for60°<a 2 ^90 c 



for 90° <a 2 < 120° 



180°-<x 



90 c 



(£ is the same for all of the slip systems within the colony.) 
30 The statistical distribution of the minimum Schmid factor can be determined 150 for the 
three slip systems as follows: 

1. Assume that the orientation of the domain is untextured 146 and thus, all orientations 
are equally likely. Therefore, C, (the angle of the slip plane normal) is uniformly 
distributed between 0° and 90°. The Monte Carlo technique can be used to choose X for 
35 each domain. 
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2. The angle between the slip direction and the loading axis Qi) is uniformly distributed 
between £ and the complement of £ as shown in Equation (27). Again, Monte Carlo can 
be used to choose X for each domain 152. 

3. Equation (30) is used to determine a t 154. 

5 4. Equation (28) is used to determine Oj and a 3 for the other two slip systems 154. 

5. Equations (31) are used to determine ^ and for the other two slip systems 154. 

6. Equation (26) is used to define M for each of three slip systems 156. 

7. The reciprocal Schmid factor for the colony is the minimum M from each of the three 
slip systems 158. 

10 8. Steps 1 through 7 are repeated many times to determine M for many colonies 160. 

^ is interesting to note that X, is always greater than X, or J^, thus ^ will never produce the 

minimum M and does not need to be considered in the above analysis.) 

The prediction from the above analysis 162 is shown in Figure 25 for 100,000 

repetitions of steps 1 through 7. The reciprocal Schmid factor for fully lamellar TiAl has a 
15 minimum value of 2 which occurs when X and £ are 45°. The maximum value is infinity. 

This is because M is high for £ near 90°. The mean value is about 3.4. 

It is to be understood that while illustrative embodiments of the invention have been 

shown and described herein, various changes and adaptions in accordance with the teachings 

of the invention will be apparent to those of ordinary skill in the art. Such changes and 
20 adaptions nevertheless are included within the spirit and scope of the invention as defined in 

the following claims. 



25 



30 



35 
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